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ABSTRACT 

Context. Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, 
and water, H2O. These high abundances would imply large cooling rates, leading to relatively short time scales for the evolution of gravitationally 
unstable dense cores, forming stars and planets. Contrary to expectation, the dedicated space missions SWAS and Odin found typically only very 
small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. 

Aims. Only toward p Oph A did Odin detect a very weak line of Ot at 1 19 GHz in a beam size of 10 arcmin. Line emission of related molecules 
changes on angular scales of the order of some tens of arcseconds, requiring a larger telescope aperture such as that of the Herschel Space 
Observatory to resolve the O2 emission and to pinpoint its origin. 

Methods. We use the Heterodyne Instrument for the Far Infrared (HIFI) aboard Herschel to obtain high resolution O2 spectra toward selected 
positions in the p Oph A core. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical 
cloud models. 

Results. The Nj = 33 - I2 line at 487.2GHz was clearly detected toward all three observed positions in the pOphAcore. In addition, an 
oversampled map of the 54 - 34 transition at 773.8 GHz revealed the detection of the line in only half of the observed area. Based on their 
ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (i 50 K) adjacent to a much colder region, where 
temperatures are below 30 K. 

Conclusions. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the 
temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on 
the derived gas temperatures, but fall into the range N(0 2 ) = 3 to ; 6 x 10 15 cm 2 . Beam averaged O2 abundances are about 5 x 10~ 8 relative 
to H2. Combining the HIFI data with earlier Odin observations yields a source size at 1 19 GHz of about 4 to 5 arcmin, encompassing the entire 
p Oph A core. We speculate that the general very low detection rate of 2 sources might be explained by the transient nature of interstellar O2 
molecules. 

Key words. ISM: abundances - ISM: molecules - ISM: lines and bands - ISM: clouds - ISM: individual objects: p Oph A SM 1 - Stars: formation 



1. Introduction 

Despite the universal importance of oxygen, its molecular form, 

02, is an elusive species of the interstellar medium (ISM). This 
molecule was thought to be one of the main regulators of the en- 



ergy b alance in the ISM, as summarised by |Goldsmith & Langer 
(1978). Consequently, large efforts, both from ground (0 15 0) 



and space, had been devoted to obtain quantitative estimates of 
its abundance. The historical account of this essentially fruitless 
"02-struggle" was recently review ed by|Gold smith et aL~|(|201 1). 
Prior to He rschel, bot h SWAS ([Melnick et al.|2000[|Goldsmith| 
et al.|2000) > and Odin ( [Nordh et. al|2003||Pagani et al.|2003) had 
already shown that abundances for both O2 and H2O assumed by 
Goldsmith & Langer| ( |1978| l were much larger than actual val- 
ues in the ISM and, with the exception of the p Oph cloud, none 
of the O2 lines were detected anywhere in the ISM. Goldsmith 



unusual velocity and was atypically broad. Pagani et al. (2003) 
showed that, on the basis of more sensitive O2 1 19 GHz obser- 
vations with Odin, this was an erroneous result. This line was 
finally detected by Odin, at the correct LSR-velocity and with a 
plausible, narrow line width ( [Larsson et al.|2007| ). Here we re- 
port the detection of two more O2 transitions in pOph A, viz. at 
487 GHz and at 774 GHz, respectively. 



Besides being relatively nearby ( 120-130 pc, Lombardi et al.| 
2008; |Snow etal.|2008|rMamajek|2008||Loinard et al.|2008| ) the 
p Oph cloud distinguishes itself from other low-mass star form- 
ing regions in that it exhibits evidence (e.g., in C ls O line emis- 
sion) of gas at relatively high temperatures (T > 20 K) over ex- 
tended regions with high column densities (7Y(H2) » 10 22 cm~ 2 , 



et al. ( 2002| ) announced the tentative detection by SWAS of the 
O2 487 GHz fine in p Oph A. The claimed signal appeared at an 



Liseau et al. 2010|). In addition, pOphA displays an interest- 



* Based on observations with Herschel-YHFl. Herschel is an ESA 
space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from 
NASA. 



ing chemistry: For example, doubly deuterated formaldehyde 
(D2CO) and hydrogen peroxide (H2O2) molecules, rarely seen 
elsewhere in the ISM have been found here (Bergman et al. 
2011a, b). The overall impression is that the observable abun- 
dance of many species is the result of surface reactions on dust 
grains, a process which may also pertain to the production of 
oxygen molecules. 
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Fig.l. Herschel beams of 2 487 GHz (44") superposed onto a grey scale C 18 0(3-2) map of pOphA ( Liseau et al.|2010 i. The + 
and x symbols designate the positions of the H- and V-polarization beam centers of HIFI, respectively, with a separation of 6'.'6. The 
heavy circles show the positions of the sources Ol, 03 and 04 (see Table[2]), with the corresponding off-source reference positions 
shown by the displaced symbols to the east and west, respectively, and identified by their observation numbers. The origin is at 
16 h 26 m 24 s 6, -24°23'54" J(2000.0) and offsets in R.A. and Dec. are in arcsec. At the distance of 120pc, 40" corresponds to about 
0.02 pc. The star symbols identify the positions of the star S 1, the young stellar object GSS 30 and the Class source VLA 1623, 
where the contours depict red- and blueshifted CO (3-2) emission of its bipolar outflow. The filled squares identify the known 
Herbig-Haro objects HI through 5 and the direction to the dominating B-star HD 147889, about 15' away, is indicated by the heavy 
arrow. 



Table 1. O2 transitions and instrumental characteristics. 



Transition" Frequency* Wavelength Energy Observatory/ Beam width Efficiency Reference 





- jvjr 


y (MHz) 


A (//in) 


E up /k (K) 


Instrument 


HPBW (') 


TJmb 


h 


- 1 


118750.341 


2526 


5.7 


Odin 


9.96 


0.91 


h 


- h 


487249.264 


616 


26.4 


SWAS 


3.5 x 5.0 


0.90 












HIFI 


0.73 


0.757 


5 4 


-3 4 


773839.512 


388 


60.7 


HIFI 


0.47 


0.753 



|Frisk et al.|(2 003 1 
|Melmck eta^O OOO 
iKoeltsema et at 



300} 
(2012) 



Koeltsema et al. 



2012) 



8 An energy level diagram is provided by, e.g., |Goldsmith et aL] ( |201l) . 

* Frequencies are from the JPL-catalogue http: //spec . jpl .nasa. gov/cgi-bin/catform. 



The motivation for the Herschel observations was, firstly, 
to pin down the precise location of the O2 source inside the 
large, ten-arcminute beam of Odin. Secondly, to add observa- 
tions of the O z 487.25 GHz (3 3 -l 2 ) and 773.84 GHz (5 4 -3 4 ) 
lines, which have upper level energies E up /k — 26 K and 61 K, 
respectively, to the Odin data of the 118.75 GHz (li-lo) tran- 
sition (E up /k = 6K). These observations should enable us to 
learn about the nature of the O2 source in the p Oph cloud. This 
would then be compared to the physical and chemical condi- 
tions of other locations in the general ISM and potentially iden- 
tify the characteristics and time scales of regions containing O2 
molecules. 



The paper is organised as follows: in Sect. 2, our Herschel- 
HIFI observations and their reduction are discussed in consid- 
erable detail and our results are presented in Sect. 3. These re- 
sults are discussed in Sect. 4, where the spectral line formation 
is analysed under different assumptions. Finally, in Sect. 5, our 
main conclusions are briefly summarised. 



2. Observations and data reduction 



Herschel is a space platform for far infrared and sub-millimetre 
observations (Pilb ratt et al.||2010 i. It is orbiting the Sun about 
1.5 million kilometres beyond the Earth (at L2) and its 3.5 m 
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Fig. 2. Averaged Double Side Band (DSB) spectra of H- and V-polarizations for the three positions Ol, 03 and 04 in the pOph A 
cloud core, shown for the range 484.65 GHz to 489.85 GHz. The spectral lines are found exclusively in the Lower Side Band (LSB) 
and the LSR velocity, i/lsr, along the abcissa is in km s and the antenna temperature Tp, along the ordinate in mK. a: 5 GHz wide 
band of the WBS around the 487 O2 line (-1600, +1600 kms~'). b: The WBS spectrum at expanded scale, showing the O2 line 
at 487 GHz, which is detected in all three positions. Inspection of the individual on-off pairs suggests that the apparent absorption 
features next to the lines are not real but simply due to the noise, c: Similar to b but for the corresponding HRS data. 

Table 2. Source designations and coordinates for p Oph A, with positions observed in O2 in bold face. 



RA (h m s) J2000.0 Dec (° ' ") J2000.0 Source Designation 



Reference 



16 25 24.32 


-24 27 56.57 


HD 147889 


16 26 17.5 


-24 23 13 


H4, HH313B 


16 26 19.0 


-24 23 08 


H5, HH313A 


16 26 21.36 


-24 23 06.4 


GSS 30, El 21 


16 26 25.7 


-24 23 24 


Ol 


16 26 25.7 


-24 23 57 


02 


16 26 26.1 


-24 23 14 


N 2 H + Nl b (FWHM=0.29kms-') 


16 26 26.38 


-24 24 31.0 


VLA 1623 


16 26 27.1 


-24 23 30 


N5, 850 jum 


16 26 27.2 


-24 24 04 


Deuterium peak 


16 26 27.3 


-24 23 28 


SM1N, 1.3 mm 


16 26 27.9 


-24 23 26 


03, P2, C 18 (3 - 2) 


16 26 27.9 


-24 23 57 


04, P3, C 18 (3 - 2); SM 1, 1.3 mm 


16 26 34.19 


-24 23 28.2 


S 1,GSS 35, El 25 



SIMBAD: http://simbad.u-strasbg.fr/simbad/ 



Dent et al. 


( 1995 1, (Jaratti (jaratti et al. (2U06) 


Dent et al. 


Jl995l,KJaratti Uaratti et al. pOOo) 


Urasdalen et al. 


j!973|, 


Llias 


1 1978 



this paper 



di Francesc o 
Andre et al 



:o et al. lt 
Jdl990l 



(2004) 



di hrance'sco et al. (2004) 

Bergm an et al. ( 201 lab 
■ ■ . i~i 1 1 1 1 J 1 1 1 1 * 



Johnstone et al. 



(2000) 



Motte et al. 
this paper, 



this paper, 



1998 



Liseau etal.|(2010) 



Grasdalen et al 



Liseau et al. (2010), Motte et al. ( 1998) 
]fT973 1, Elias ( 1978 ~ 



primary dish is radiatively cooled to its operational equilibrium 
temperature of about 85 K. The scientific instruments are placed 
in a liquid helium filled cryostat, limiting its cold lifetime to 
roughly 3.5 years. One of the three onboard instruments is the 
Heterodyne Instrument for the Far Infrared (HIFI, |De Graauw| 
et al. |2010 1 with continuous frequency coverage from 480 to 
1250 GHz (/LI 624 - 234 pm) in 5 bands. In addition, the fre- 
quency range 1410 to 1910GHz (/L1213 - 157pm) is covered 
by bands 6 and 7. The spectral resolving capability is the highest 
available on Herschel, i.e. up to 10 7 corresponding to 30 ms -1 , 
required to resolve the profiles of very narrow lines of cold 
sources, having temperatures of the order of 10 K or less. 



Spectral line data and relevant characteristics of the instru- 
ments discussed in this paper have been compiled in Table[T] 



2.1. 487 GHz observations 

The 487 GHz observations with HIFI were performed on oper- 
ating day OD673 (2011 March 18). The three positions Ol, 03 
and 04 at the centre of the p Oph A core (Table[2]i were observed 
in double beam-switch mode (DBS), with beam throws of 3' 
in roughly the west and east directions (with a position angle 
of 101° E of N) and identified in Fig.[T]by the numbers 6031 
and 6032, respectively. For each position, the relative observing 
time spent on the east- and westward on-off pairs was 7.3 hr, and 
the total programme execution time was 21.9 hr. The DBS mode 
generally produces the highest quality data with HIFI and since 
the O2 signal is known to be weak, this was the observing mode 
chosen. Switching by only three arcminutes inside a molecular 
cloud is generally considered not a sensible option, as this po- 
tentially results in cancellation of the source signal. However, 



3 



R. Liseau et al.: Multi-line 




Fig. 3. Outline of the oversampled 6x6 regular map with 10" 
spacings of the O2 774 GHz observations. At this frequency, the 
beam of the 3.5 m Herschel telescope has a size of 28" (small 
blue circles). The beam centers for the H- and V-polarisations 
are separated by 4'.'5. Larger (black) circles correspond to the 
487 GHz beam positions (44"; cf. Fig[T]i and as before, North 
is up and East is to the left. Background image is courtesy of 
Sp/fzer-IRAC-MIPS (blue: 3.6 //m, green: 8/im, red: 24 /mi). 



of 2 toward p Oph A 



Table 3. O2 line observations of p Oph A. 



Position/ 


mb 







FWHM" 


10 3 (T mb dv 


Average 


(mK) 


(kms- 1 ) 


(km; 




(Kkm 






487 GHz 


Nj = 


3 3 -l 2 










01 HRS 


14.4 ± 


1.8 


3.05 


±0.10 


1.63 ± 


0.23 


19.8 ± 


2.4 


01A b 


24.2 ± 


3.3 


2.79 


±0.05 


0.46 ± 


0.08 


15.5 ± 


1.7 


r\ 1 o 


19.6 ± 


3.5 


3.58 


±0.04 


0.41 ± 


0.09 


8.7 ± 


1.7 


03 HRS 


7.6 ± 


2.2 


3.59 


±0.12 


0.80 ± 


0.27 


5.7 ± 


1.3 


04 HRS 


11.5 ± 


2.5 


3.75 


±0.07 


0.84 ± 


0.16 


13.0 ± 


1.6 


Ol WBS 


9.4 ± 


2.0 


2.89 


±0.20 


1.91 ± 


0.46 


18.4 ± 


1.6 


oi_i 


9.6 + 


1.4 


2.70 fix 


0.90 


fix 


9.2 ± 


1.4 


01_2 


6.9 + 


1.4 


3.60 fix 


0.90 


fix 


6.6 ± 


1.4 


03 WBS 


5.7 ± 


3.4 


3.66 


±0.26 


0.87 ± 


0.60 


5.3 ± 


1.5 




10.3 ± 


2.6 


3.78 


±0.12 


0.99 ± 


0.28 


10.3 ± 


1.5 


<03, 04>WBS 


7.8 ± 


3.0 


3.69 


±0.20 


0.90 ± 


0.47 


7.8 ± 


1.5 




774 GHz 


Nj = 


54-34 










<01, 02>HRS 


22.3 ± 


9.9 


3.59 


±0.03 


0.57 ± 


0.10 


13.2 ± 


4.0 


Ol WBS 


17.9 ± 


4.2 


3.66 


±0.07 


0.61 ± 


0.16 


11.4 ± 


2.9 


r\i \X7R c 
UZ WJdo 


14.5 ± 


3.8 


3.59 


±0.09 


0.72 ± 


0.22 


10.9 ± 


2.8 


<01, 02>WBS 


15.9 ± 


2.8 


3.62 


±0.06 


0.67 ± 


0.17 


11.2 ± 


2.0 


03 WBS 














3.6 ± 


2.9 


04 WBS 














-1.6± 


3.0 


<03, 04>WBS 














0.8 ± 


2.0 




119 GHz 


Nj = 


li-lo 










Odin 1 


19.3 ± 


0.6 


3.5 


±0.5 


1.5 ± 


0.5 


26.7 ± 


4.4 



8 Gaussian fits: peak intensity, central velocity and full width half 

maximum of line. 
* 01_1 and 01J2 refer to two velocity components (see Sect. 3.1). 
c 10' beam centered on p Oph A ( Lars son et al.|2007) . 



a relatively large velocity gradient (>0.3kms '/arcrnrn) and/or 
two distinctly different velocity systems (Ay ~ lkms" 1 ) are 



known to exist in pOphA (e.g., Bergman et al. 201 la 1. In ad- 
dition, the O2 line is known to be narrow (FWHM ~ 1 kms 1 ), 
which taken together largely diminish the risk of signal cancel- 
lation. In fact, splitting the entire data set into two and reducing 
these two halves with only one off-spectrum used resulted in es- 
sentially the same spectral data. These off-beams were pointing 
at very different regions in p Oph A (Fig.[T]i, but no traceable off- 
beam contamination was introduced by the DBS observations. 

The receivers are double side band and to counteract confu- 
sion in frequency space, we had selected for each observation 
8 different LO-tunings inside HIFI-band la, spaced by 170 MHz 
to 260 MHz. During data reduction, this allowed sideband de- 
convolution and the identification of any "false" signal due to a 
strong feature in the other sideband, which might otherwise be 
folded over onto the weak O2 feature. 

Both the Wide Band Spectrometer (WBS) and the High 
Resolution Spectrometer (HRS) with respective resolutions 
1.10MHz and 250kHz (0.68 and 0.15 km s 1 ) were used. The 
positions of spikes and spectrometer artifacts are known and 
flagged and none of these interfered with the O2 observations. 
Initial data reduction and calibrations were performed with the 
standard Herschel pipeline HIPE version 7. 1 and the subsequent 
data analysis was made with a variety of other software pack- 
ages. The quality of the HIFI data is excellent and it was suffi- 
cient to subtract a low order baseline for each LO setting prior to 
averaging. The intensity calibration accuracy is estimated to be 
10%. 



The frequency dependent beam widths (HPBW) and main 
beam efficiencies (rj m b) were adopted from Roelfs ema et al.| 
(2012 1. Before averaging the data for a given position, the H- and 



V-polarization spectra were all inspected individually in order to 
avoid the suppression of weak spectral lines or the generation of 
artificial features. 



2.2. 774 GHz observations 

The 774 GHz observations with HIFI Band 2 were made on 
OD 583 (201 1 September 15). A 6 x 6 map, aligned with the 
equatorial celestial coordinates (Fig. [3]), was obtained with 10" 
regular spacing, oversampling the 28" beam by nearly a factor of 
three and allowing for potentially increased spatial information 
of the emission regions. 

Both the WBS (0.43 km s 1 ) and the HRS (0.097 km s" 1 ) 
were used. The 774 GHz data were also obtained in DBS mode, 
with the offset positions displaced by 3' and in nearly the same 
chopping direction as before (along pa = 98°). With 500 s on- 
source integrations at each position, the total observing time was 
12hr. 



3. Results 

In the figures showing HIFI data, brightness is given as antenna 
temperture, T A . When more than one telescope had been used, 
the data are presented on the main beam brightness temperature 
scale, r m b, where the main beam efficiency is rj m \, = 0.75 for 
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both HIFI trans itions ([R oelfsem a et al.||2012|l and 0.9 fo r the 
Odin transition ( Frisk et al.||2003| Hjalmarson et al.||2003] > and 
where we use the relation 7A/?7mb = 7mb- For completeness, we 
also provide in Table[3] together with the HIFI data, the results 
of the Odin observations. 



3.1. 487 GHz 

The reduced WBS spectra for the three positions 01, 03 and 04 
are shown for the frequency range 484.65 GHz to 489.85 GHz 
in Fig. [2^. The Upper Side Band (USB) is completely line-free, 
i.e. all line detections refer to the LSB only. The presented data 
are DSB in order to suppress the noise. Centered on the O2 
line, radial velocities span -1600 to +1600 kms~'. In addition 
to the detected O2 487.25 GHz lines near an LSR velocity of 
3kms _1 in all three positions, the spectra show a number of 
other emission features which correspond to identified molec- 
ular lines, due to, e.g., CS, CH3OH, NH 2 D and SO, to mention 
only the strongest ones. Also these spectra reveal variations in 
intensity and line width on angular scales of only some tens of 
arcsec. This remarkable behaviour of different molecular species 
toward the p Oph A core has been known for some time, and has 
recently been documented by Berg man et al.| ( [201 la) . Our multi- 
line HIFI data will not be discussed further here but be presented 
elsewhere (Larsson et al., Pagani et al., in preparation). 

Figure^ shows the WBS spectra zoomed in on the 
O2 487 GHz transition. This line too displays varying intensity 
levels in the three positions. In particular, the line toward the Ol 
position appears, by comparison, rather strong and broader than 
those toward 03 and 04. The HRS data shown in Fig. [2}; sug- 
gest the line at Ol to be a composite of two different, narrow 
components, blending into a single broad feature at the lower 
resolution of the WBS. Henceforth we will refer to this spec- 
tral component at l<lsr — f-3.6kms _1 as 01_2, whereas the one 
closer to 2.8 kms 1 is referred to as 01_1. 

We summarise the spectral characteristics of the O2 487 GHz 
line in Table[3] Within the errors, the WBS line centroids and 
widths are in agreement with the Gaussian fittings of the HRS 
data. Also, the average intensities for positions Ol J2, 03 and 04 
with the WBS (8.3 + 1 .8 mK km s 1 ) and the HRS (9.1 ± 1.9 mK 
kms~'), respectively, are consistent with each other. 

In addition, the O2 line centroid and width also confor m with 



AR.A. ["] 



values of other optically thin species, e.g. C O (3 - 2) (Liseau 



etaLl|20T0]> N 2 H + (1-0) (di Francesco et"alj|2004l ), deuterated 
molecules (Bergman et al. 2011a) and the recently discovered 
H 2 2 ( I IOOI I. Bergman cl al.|2()1 1 h). Regarding the 2.8 km s 1 
feature, we can at present not entirely exclude the possibility 
that it is due to an unidentified species. If so, the line transition 
frequency should be v = (487250.548 ± 0.098) MHz (lcr). 



3.2. 774 GHz 

The convolved and averaged 774 GHz spectra are compared to 
487 GHz spectra in Fig.^and the 6 x 6 map of the 774 GHz 
observations, convolved to the beam at 487 GHz, is displayed 
in Fig. [5] The line was not detected toward the east side of 
the cloud, i.e. toward 03 and 04, with lcr-upper limits of 
3mKkms _I , but clearly detected toward the west side of the 
cloud, i.e. toward Ol and 02, at S /N — 4 - 6. The derived line 
parameters for the 774 GHz spectra, convolved to the 487 GHz 
beam, are reported in TablepT We note that only the 3.6kms~' 
component is seen at the Ol position. 




-5 5 10 



-5 5 10 



[km/s] 



Fig. 4. The WBS spectra toward the positions Ol to 04 are 
shown for both the 774 GHz (upper) and 487 GHz (lower) lines 
in each frame. The averaged 774 GHz data are convolved to 
the 487 GHz beam of 44". Offsets relative to the origin 04 are 
shown along the upper and right-hand scales. 



The 2 774 GHz line falls in the lower side band (LSB: 
770.75-775.65 GHz). The entire observed spectral region is 
however totally dominated by the 13 CO (7-6) transition at 
771.2 GHz. This line exhibits significant variations on small an- 
gular scales which is clearly evident in individual pointings of 
the 774 GHz beam of 28" (FWHM). In the spectrum farthest to 
the northeast, the line reaches a peak T m \, of 25 K, whereas in the 
opposite corner in the southwest the intensity has decreased to 
less than 5 K. 

In addition, in the upper side band (USB: 783.15- 
786.815 GHz), the strongest and clearly detected line is due to 
the C n O (7-6) transition at 786.3 GHz (Larsson et al., in prepa- 
ration). In contrast to 13 CO, the line of C 17 has its maximum 
of almost 3 K at the Ol position. 

4. Discussion 

4.1. Observed quantities: Temperature, line width, column 
density and size of the O2 emission region 

The O2 487 GHz and 774 GHz lines are clearly detected toward 
multiple positions. From the observation of more than one tran- 
sition, some basic physical parameters of the source can be de- 
rived, where we assume that all transitions trace the same gas 
along the line of sight. In the case of the observed O2 emission, 
this should be straightforward, as detailed computations of sta- 
tistical equilibrium and radiative transfer for a multi-level system 
demonstrate this emission to be optically thin and to originate in 
conditions of local thermodynamic equilibrium (LTE). 

4.1.1. Temperature 

At Ol, the observed source extent of the 487 GHz and 774 GHz 
lines seems similar so that the ratio of the beam fillings 
fAKilfnA is probably not far from unity. In addition, the main 
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Table 4. Observed O2 line intensity ratios. 



Spectrum/ 




T, ■ 

1 kin 




Position 


7(774) 


(K) 


(cm- 2 ) 


01_1 (2.8 kms -1 ) 


> 0.9 (3o-) 


< 40 


> 2 x 10 15 


01_2(3.6kms-') 


0.58(±0.19) 


78 


5.4 x 10 15 




(+lo-) 0.77 


47 


5.0 x 10 15 




(-Icr) 0.39 


» 100 




03 


> 0.61 (3<r) 


< 70 


> 3 x 10 15 


04 


> 1.14(3cr) 


< 31 


> 6 x 10 15 


<01,03, 04> 


1.43(±0.94) c 


25 


3.4 x 10 15 




(+lo-) 2.37 


18 


4.6 x 10 15 




(-lcr) 0.49 


> 100 





" Errors are ±(hlh) V(<5/i//i) 2 + (ShJh?- 
*> For HIFI, 17^(487) = r/ mb (737) = 0.75. 

c Similar to 03 for lcr limit on 774 GHz intensity (1.47 ± 1.26). 



beam efficiencies at these frequencies are also about equal 
(Table[T|. Therefore, the observed ratio of integrated intensities 
7(487)//(774) bs is directly comparable to the theoretical one, 
7(487)//(774)th eo , obtained from LTE calculations for optically 
thin O2 radiation and which as such depends only on the gas 
temperature (Fig.|7|. 

From Table|4] we see that the measured intensity ratio for the 
487 GHz and 774 GHz lines toward 1 implies a temperature of 
nearly 80 K, with the large error allowing a large range, from 
about 50 K to very much higher than 100 K (see also Fig.|7]i. If 
the beam filling at 774 GHz were smaller than that at 487 GHz, 
the actual temperatures would be even higher. Our value can 
be compared to earlier estimates, obtained at comparable spatial 
resolution, of both the dust and the gas temperature toward Ol 
From FIR continuum measurements with a 40" beam, Harvey et 
al. ( 1979 1 obtained the dust temperature 7a = 41 K. Also with a 
40" beam, Zeng et al. ( 1984 ) derived from NH3 observations the 
temperature of the gas, T g ~ 45 K. 

For the dense spot 04 (P 3, SM 1), the 3rr upper limit on the 
774 GHz line yields a temperature strictly lower than 3 1 K, which 



would be in agreement with T g = 22 ± 3 K obtained by [Bergman 
et al. ( 201 lb| l. For the dense core p Oph A, the average tor <OI, 
03, 04> results in a value of about 20 to 25 K, which is entirely 
consistent with earlier determinations, based on a variety of tech- 
niques. Temperatures between 9 K and 49 K have been derived, 
but with most values cluste ring around 20-30 K (e.g.,|Harvey et| 
al.|1979||Lo"ren et al.|1980[|Ward-Thompson et al.|1989||Motte 
et all 19981 |Stamatellos et al.|2007||Bergman et al'|2011a||Ade| 
et al. [2011) . 



As is also evident from Fig. [5] a temperature gradient or 
two distinctly different temperature regimes, i.e. with larger than 
about 50 K and strictly lower than 30 K, exist within the bound- 
aries of our limited map of roughly one square arcminute (see 
also Figs.[TJand[3]l. 



just barely resolved with the WBS (5v = 0.43 km s _1 ). However, 
the average 774 GHz spectrum <01, 02> is well resolved with 
the HRS and the line has a comparable width to what is observed 
with the WBS. What is remarkable is that the lines in the cold 
regions 03 and 04 are observed to be broader than those in the 
warmer Ol and 02. 

The non-thermal contribution to the observed width of an 
unblended line is u non _ th = (FWHM 2 - v^) 1 ' 2 , where all radial 
velocities are expressed in terms of their full width half max- 
imum value. For T in the range 18 K to 30 K, 1% = 0.27 to 
0.35 km s -1 , yielding c n on-th = 0.72 to 0.80 km s -1 for the cold 
region including 03 and 04 (Table|4|i. Hence, we find that there 
the line broadening is totally dominated by the non-thermal ran- 
dom motions, being more than twice as high as the thermal ones. 
On the other hand, in the warm region Ol and 02, the thermal 
motions dominate over the non-thermal ones (0.57 km s _1 versus 
< 0.44 km s~' for T = 78 K) or being comparable at the lower 
limit temperature of 47 K. This could suggest that, if the non- 
thermal motions are identified with turbulence at Ol and 02, 
much of this turbulence has been dissipated and converted into 
heat. 

The 02 source is seen in projection close to the outflow from 
VLA 1623 and judging from its position (Figs.[T]and[3]), it is pos- 
sible that the gas at 02 is enriched by shock processed mate- 
rial. At the border of the outflow, a broader line might be ex- 
pected due to enhanced turbulent motions. In fact, the 774 GHz 
line width appears (marginally) broader than that toward Ol 
(Table|3]l, but the uncertainties are large and the significance is 
low. 



4.1.3. Source size 

Having determined the temperature, we can use the theoretical 
line ratio for another transition and equate it to the ratio of the re- 
spective beam fillings, e.g. R = (/ll9//487)theo(r)/(/ll9// 4 87)ob S = 

/487//119. The relevant 119 GHz data have been adopted from 
Larsson et al. (2007 1 and are compiled in Table[3] and all 
lines are displayed on the r m b scale in Fig. [6] For a Gaussian 
beam at 119 GHz, the beam filling factor is given by fu9 ~ 

^119/(^119 + ^odin-*' wnere 0119 is me size of the source and 80dm 
that of the beam at 119 GHz. Re-arranging, this reads #119 w 

#odin/ VW/119 - 1 = 0<W If 4ii - 1, where 6> din = 10'. 

The 487 GHz source seems to be extended on the scale of at 
least the beam width, i.e./^ > 0.5, and for the average values of 
<01, 03, 04 > in Table|4] we find a size of the Odin source of 2 4 
to 5 arcmin, for gas at 18 to 25 K, i.e. for (/ii9//487)theo(T) =11 
to 8, respectively. The solution for a much warmer gas (127 K 
for "-lcr" in Table[4]i ought to be excluded, as the calculated 
source size would exceed 9', whereas the high temperature trac- 
ing 774 GHz line has been detected only over a limit ed region 
(< 1'). On the basis of CS (2-1) mapping observations, Liseau et 
. ( 1995 ) determined the deconvolved angular size of the dens- 
est regions of the p Oph A core as FWHM ~ 4' («h ~ 10 5 crrT 3 ), 
a p Oph A size consistent with the HIFI and Odin observations 



of O2. It appears that the entire p Oph A core is faintly glowing 
in the 119 GHz line. 



4.1 .2. Line width 

Comparing the observed line widths, FWHM, of Table[3] with 
the resolutions offered by the HRS and WBS, we see that the 
HRS data for the 487 GHz line are all well resolved (Sv = 
0.15 km s 1 ), whereas the 774 GHz line toward Ol and 02 is 



4.1 .4. Column density 

For the warm region Ol, an O2 column density of 5 x 10 15 cm -2 
is derived, whereas for the cold 04, N(Oo)> 6 x 10 15 cm~ 2 
(TableQ. These estimations are based on the assumption of 
optically thin emission in LTE at the kinetic gas temperature 
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Fig. 5. Left: The regridded 6x6 map of the WBS 774 GHz spectra, weighted by the Gaussian profile of the 487 GHz beam of 44". 
The darker area shows roughly the region inside that Half Power Beam Width (HPBW). The vlsr and 7a scales are indicated in the 
upper right corner. Right: Same as to the left but shown as an image of the spatial distribution of the integrated intensity, jT^dv, 
with observed positions shown as plusses. The circles show the 487 GHz and 774 GHz beams, respectively. 



7\in, i.e. for unit beam filling the O2 column is given by 



-2 



N(Q 2 ) = fT mh dvx <D(r kin ) cm 



when the integral is expressed in K cm s 1 and where 



1 tr 



gvfflda) 



Fv(7- ki „)-F„(r bg ) 



Q(T kin ) exp(r u p/rki n ), 



in K -1 crrT 3 s (cgs units throughout). Here, the transition tem- 
perature is T tr = hv/k, the quasi-Planck function is F V (T) = 
T tr /(e Tl 'l T - 1), the partition function is Q(T), the upper level 
energy is 7\, p = E up /k and the other symbols have their usual 
meaning. The background temperature is that of the cosmic mi- 
crowave background, i.e. 7bg = TcMBnThe frequencies were 
adop ted from |Drouin et aL] ( |2010| l ana the Einstein A-values 
from |Marechal et al.| ( |199"7] l (see also |Black & Smith|1984) . 

If on the other hand <01, 03, 04> can be taken as a rep- 
resentative for p Oph A as a whole, column densities lower than 
5 x 10 15 cirT 2 would be indicated. Although there is a consider- 
able range in gas kinetic temperature, the O2 column densities 
seem not to vary much. However, with the observed inhomo- 
geneity of the emission in mind, the sigificance of such an aver- 
age could be questioned. 

4.2. 2 abundance 

Once the O2 column density is known, the line of sight average 
O2 fractional abundance, with respect to H2 or H-nuclei, can be 
measured if the column density of H2 molecules is available. The 
hydrogen column density along the appropriate lateral and depth 



1 In our numerical radiative transfer calculations (Sect. 4.1), also the 
dust continuum is included. These calculations make use of the collision 
dataof |Lique|p010 i. 



scales is not well known from observation. However, we can get 
an approximate idea from the C 18 observations by 



al. 



Liseau et 



( 201 0| ), where it was suggested on the basis of the line pro- 
files that the 7=3-2 line was a reasonably good proxy for the 
2 1 19 GHz line observed with Odin ( |Larsson et aT]|2()07l ). If 
the O2 molecules are indeed cospatial with those of CO and its 
isotopic variants, their optically thin lines could be used to infer 
the H2 column densities relevant also for O2. We derive "local" 
C 18 column densities taking C I8 intensity and O2 tempera- 
ture variations into account. The gas and dust peaks coincide 
spatially (see also Table[2]i, i.e. there is no indication of signif- 
icant CO fre eze out (at levels high er than a factor of 2 to 3, 
see Fig. 10 of Bergman et al.|201 la 1. The C ls O column density 
varies by less than 20% for 20 K < T^ n < 50 K. Assu ming that 
C 18 0/H 2 = 1.5 X IP" 7 (likely wit hin about 30%, e.g. |Wannier 
|et al.|1976||Goldsmifh et al.|2000| i, H 2 columns are found to be 
1 x 10" ciir 7- 



Hence, for the O2 column densities in Table|4] we find that 
X(C>2) = 5 x 10~ 8 relative to H2 in the warm gas toward the west, 
at 01 (and likely also toward 02), and higher than that toward 
the east, at the colder 04. The C ls O lines are slightly optically 
thick, implying that these H2 column densities are lower limits 
and need to be adjusted upwards by a factor of the order of t/( 1 - 
e~ T ) i 2, for r ~ 2. Hence, these O2 abundances are upper limits, 
i.e. the fractional abundance of O2 in the gas phase is less than 
indicated above. 



A way to obtain a more detailed picture of the source is to 
use detailed theoretical models to predict the O2 abundance or 
column and to compare these models with the observations. This 
will be the topic of the next section. 
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Fig. 6. Averaged data for the HIFI O2 observations (774 and 
487 GHz) and the Odin 119 GHz spectrum toward pOphA. 
For proper comparison, these spectra are given in the T m \, 
scale and for reference, a vertical dotted line is shown at 
vlsr = +3.5kms~ 1 . 



4.3. Chemical models 



Specifically for Q 2 (and H 2 Q), Holle nbach et al.| ( |2009] l 



constructed theoretical models of Photon-Dorninated Regions 
(PDRs, Hollenbach & Tielens 1997 1, invoking also grain surface 
chemistry, including freeze-out and desorption of the species. 
These new models are models of the whole cloud, to arbitrary 
values of the extinction, and showing that most of the gas phase 
O2, OH and H2O is found at values of Ay ~ 3-10 mag, where the 
details depend, as in a conventional PDR, on the ratio Gq/«h- In 
many aspects, these models are similar to the traditional dense 
PDR models including chemistry and gas heating which have 
been developed to explain the abundances and column densities 
of other molecules (including radicals and ions) at intermedi- 
ate depths into the clouds (e.g ., Fuente et al.|1993 Sternberg & 



Dalgarno 1995) |Jansen et al.||1995[ |Simon etal 



1997 Pety et 



al. 2005|l. The current models are extended even deeper into the 
cloud where external UV photons play no role anymore in the 
chemistry, only cosmic -ray induced photons. 

In the direction of pOphA, a number of spectroscopic in- 
dicators suggest the presence of a PDRp] |Liseau et al.| ( |l999| l 



2 Atomic recombination lines of carbon, 12 C, and sulfur, 32 S, 
(Cesarsky et al. 1976; Pankonin & Walmsley 1978), fine structure lines: 



[C I] 492 . 2 GHz, (TKamegai et al.||2003| >, [CI] 809.3 GHz JKamegaTet 
[aLl [2003] |Kulesa et al.||2005| >, [C II] 1 57 pm \Yui et al.|[l993| |Liseau 
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Fig. 7. Upper: LTE computations of the O2 column density 
N(02) as function of the temperature T for an extended homo- 
geneous source, to yield a line intensity of 1 .0 Kkm s _I . For op- 
tically thin emission, intensity and column density are linearly 
related. 1 19 GHz: blue dots; 487 Ghz: black solid line; 774 GHz: 
red short dashes. Lower: The ratio of the 487 GHz integrated 
line intensity to that at 774 GHz is shown by the curve for a range 
of temperatures. For Ol (= Ol _2), the error shown by the dashed 
line formally extends to the unplausible result of nearly 10 3 K. 
For 04, the limit symbols are for the lower limit on the intensity 
ratio and the upper limit on the temperature, respectively. 



inferred that the PDR is situated on the rear side of pOphA. 



as seen from the Sun, and excited by the two stellar sources 
HP 147889 and S 1 (see Fig.|l| and Table|2|> by their combined 
UV fields, corresponding to Go ~ 10 z . For the densest parts of 



|eTaT|T999l >, [O I] 63 pm <|Ceccarelli et aLp997l |Liseau et al.|[l999> , 
[O I] 145 pm (Liseau et al.|1999||Liseau & lusttanont |2009f and mid- 
infrared PAH emission: 1 Liseau & lusttano nt |2009) . 
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the [C II] emitting p Oph A-PDR these authors derived empiri- 
cally «h = 5 x 10 4 crrT 3 and over more extended regions, 1 to 
3xl0 4 ctrT 3 . 



The model of Hollenbach et al. ( 2009 ) predicts columns of 



O2 between about lO 15 and a few times 10 cm depending 
on the grain temperature in the 02-emitting region, with higher 
columns in regions with warmer dust. Warmer (7^ > 20 K) 
dust produces more O2 column because O atoms evaporate from 
grains before they can react with atomic hydrogen sticking to 
the grains. This increases the elemental gas phase O available 
to make O2, whose abundance in the pOph model can reach lo- 
cally above 10~ 5 for = 25 K. The results are quite insensitive 
to the gas temperature in the 02-emitting region for the range 
10 K< T g < 100 K. 

However, the [Hollenbac h et al.| ( |2009[ ) models predict gas 
temperatures of only ~ 20 K, as would all thermochemical mod- 
els with only cosmic ray heating and heating by the penetrating 
FUV. Gas temperatures as high as 80 K, as derived here, sug- 
gest other heating mechanisms, such as slow shocks or turbulent 
dissipation (e.g., Golds mith et al.|2 010). Our observations show 
a slight increase in O2 column in cooler gas. If the dust is also 
cooler in the cooler gas region, as would be expected for the de- 
rived densities above 10 4 and, centrally, higher than 10 7 cm~ 3 , 
this contradicts the predictions of this particular model. 

On the other hand, geometry also plays a role in deter- 
mining the column of O2 observed along the line of sight. 
The [Hollenbach et al. (2009 ) models assume an illuminated 
slab face-on. However, dumpiness and more edge-on geometry 
could increase the observed column in some beams, and could 
explain the slightly increased column of O2 in the region with 
somewhat cooler gas. We conclude that these models are mod- 
estly consistent with the observations, but likely require some 
additional gas heating in the 02-emitting region. We defer de- 
tailed thermochemical modeling of these p Oph observations to 
a future paper. 



4.4. Source(s) of O2 or the dearth of O2 

The discussed models could be expected to apply quite generally 
to all externally illuminated molecular clouds developing PDR- 
like parts. It is surprising therefore that O2 has been detected 
only toward the dense cloud core p Oph A ( Larsson et al.| 2007 



and this work) and a spot in OMC-1 (Goldsmith et al.|201 l| l, in 
spite of numerous searches toward other regions. For instance, 
surveying twenty sources, Goldsmith et al. (2000) obtained lim- 
iting O9 abundances 7V(02)/W(H2) of a few times 10" 7 , whereas 



Pagani et al. (2003) derived limiting O2 columns of 10 15 to a 
few times 10 cm for nearly a dozen objects. These surveys 
included dark clouds, regions of low mass star formation in the 
solar neighbourhood and of high mass star formation at different 
galactocentric distances, i.e. these surveys sampled a wide range 
of cloud properties and energetics, from quiescent cold clouds 
to turbulent shock heated gas due to mass infall and molecular 
outflows. Yet, most escaped detection, although their limiting O2 
column densities fall within the range predicted by the theoreti- 
cal models. 

One possibility could be that this apparent mismatch is due 
to observational bias, such as source distance and beam dilution. 
For instance, the limit on the O2 column in Ori A determined by 
Paga ni et al.| (f2003) is entirely consistent with the HIFI beam- 
averaged detection obtained by Goldsmith et al. ( 201 lj ), if proper 
correction for dilution in the Odin-beam is applied. We recall 
that the detected spot is not toward the Orion-PDR ( Meln ick~et| 



|al.|2012| ), as one might have expected, but close to the outflow 
in the vicinity of IRc 2. 

However, the difference in distance aside, with differing 
Go I 'rin, pOphA and Ori A represent very different environ- 
ments, viz. cold core deuteration versus hot core grain chem- 
istry (Bergman et al. 2011a Bisschop et al. 2007 Ceccarelli 



et~aLl|2007| |Ratajczak et al.||2011) . Star forming regions like 
IRAS 16293 in Ophiuchus, NGC 1333 in Perseus or the Serpens 
core near S 68 are more akin to the conditions in p Oph A, yet 
no O2 detections have been reported. These complexes are at 
the same or within at most twice the distance to pOphA, so 
that beam dilution issues should be much reduced. It is conceiv- 
able that the viewing geometry of p Oph A could be particularly 
favourable for the observability of O2 and that this aspect is not 
shared by other regions. 

However, a clear difference is suggested by the fact that hy- 
drogen peroxide (H2O2) has so far only been detected in p Oph A 
(Bergman et al.|2011b| ). H2O2 is thought to be a clear indicator 
of grain surface chemistry and relatively abundant at about 30 K 
(e.g., Cazaux et al.|2010] ). The simultaneously detectable abun- 
dance of O2 and H2O2 may perhaps be used to infer the evolu- 
tionary chemical state of the cloud, hence its "age", if interstellar 
gas phase O2 is relatively transient ( Du et al.|2012 ). In that case, 
it could explain the paucity of known O2 line emitters. 



5. Conclusions 

• We have successfully detected the 33 - I2 transition of O2 at 
487 GHz with HIFI onboard Herschel toward three positions 
in the p Oph A core. In addition, we have also obtained a 6x6 
oversampled map with 10" spacings in the 54 - 34 transition 
at 774 GHz. 

• The telescope pointings were toward two cold, high den- 
sity cores, i.e. 04 [C 18 OP3 = SM 1] and 03 [C 18 OP2, near 
SM IN], and two positions 30" west of these, labelled 02 
and Ol, respectively. All three observed positions (Ol, 03 
and 04) were detected in the 487 GHz line, with significant 
variations of the line intensity on scales as small as 30". At 
774 GHz, emission was detected from only the western part 
of the map, including Ol and 02. 

• Combining the 487 GHz and 774 GHz data leads to the col- 
umn density of N{02) > 6 x 10 15 crrT 2 , at the 3cr level and at 
T < 30 K, in the high density region 03 and 04, whereas in 
the warmer region of Ol and 02, 7V(02) = 5.5 X 10 15 crrT 2 
{T > 50 K). There, standard analysis yields an abundance 
of Af(02)//V(H 2 ) ~ 5 x 10~ 8 in the warm gas and somewhat 
higher in the cold region. 

• This result is in agreement with X(02) based on the obser- 
vation of the 119 GHz line with Odin ( |Larsson et al.|2007| ), 
assuming that the O2 source size has an angular extent of 
about 5'. 

• Comparing our results with the predictions from theoretical 
chemical models, we find good agreement between observed 
and predicted O2 column densities. Less convincing agree- 
ment regards theoretical and observed temperatures, the rea- 
sons for which are not clear. 

• The question why O2 is such an elusive molecule in the ISM 
is still unanswered. In the special case of p Oph A, there is 
some evidence suggesting that detectable amounts of gas 
phase O2 might be a relatively transient phenomenon, which 
could explain why interstellar O2 is generally not detected. 
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